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Abstract 



CN ' The properties of pairing correlation in nuclear matter are investigated by using 

. various versions of Skyrme forces. Truncation of states involving pairing correlation, 

| necessary due to zero range nature of the Skyrme force, is discussed in detail. A plateau 

appears in pairing gap versus cutoff for each force. We propose to choose the cutoff 
parameter in the middle of the plateau so that the parameterization is independent of 



nuclides. 



The mean field theory with the Skyrme force [|l], [2| has been successful in describing 
properties of the atomic nucleus. Bulk features such as binding energies and radii have been 
reproduced well and, in addition, more detailed properties such as fission barriers, giant 



resonances, single-particle spectra and low-lying coherent particle-hole and particle-particle 

h : 

excitations have been studied. Various versions of the force have been constructed for specific 
purposes. 

Dobaczewski et al. || proposed the SkP force in order to describe the pairing properties. 
While pairing correlation contributes to various properties of nuclei, one of the most inter- 
esting topics is the formation of skins and halos in neutron excess nuclei. In these nuclei 
the neutron Fermi energy rises close to but below zero so that unoccupied neutron orbitals 
just above the Fermi level are unbound. Therefore the pairing interaction scatters particle 
pairs from bound states to continuum states. This process can have sizable amplitude in 
the peripheral region and the formation of neutron skin or halo can be enhanced. Pairing 
correlation in neutron excess nuclei is also discussed in other articles, e.g. ref. [|J. 
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In the Hartree-Fock (HF)+BCS method, the mean field is solved self-consistently but 
the pairing field is separately treated. The pairing interaction is usually taken to be so 
called seniority force or the monopole force. The force is independent of the orbital so that 
gap equations are decoupled from the HF field and the whole scheme becomes extremely 
simple. For each nucleus, the strength of matrix element G is determined by, e.g. fitting to 
experimental odd-even mass staggering. The matrix element G must, however, be dependent 
on quantities other than N and Z, e.g., the density and deformation. In order to deal with 
the pairing gap more properly, it is desirable to utilize interactions independent of nuclides 
like the Skyrme force as a more fundamental approach. 

The treatment of HF+BCS including continuum states may suffer from unlocalized den- 
sity due to unphysical "particle gas" surrounding the nucleus. As shown by Dobaczewski 
et al. P|, Hartree-Fock-Bogoliubov (HFB) theory is free from this difficulty In this theory 
both the mean field and the pairing field are treated on an equal footing and determined 
simultaneously from the same interaction. 

Though HFB has been performed with finite range Gogny force 0, it is practical to 
employ HFB with the Skyrme force in describing neutron excess nuclei. In using Gogny 
force the expansion in terms of the harmonic oscillator basis provides a poor approximation 
to continuum states, even when a large number of oscillator basis is taken into account. 
Thus coordinate-space treatment is necessary for neutron-excess nuclei. The Skyrme force 
has the possibility of working in coordinate space. This is the reason why the Skyrme force 
is employed as a pairing force in describing nuclei far from stability ||. 

The consistent usage of the Skyrme force is also preferable in the application of Generator 
Coordinate Method [|]. 

In the present work we investigate the properties of pairing correlation of Skyrme forces. 
Firstly, the formalism and notation are presented. Secondly, the strength of pairing correla- 
tion in nuclear matter is calculated by making use of various versions of Skyrme forces. The 
main subject of this paper is the dependence of the pairing gap on the cutoff parameter, 
which is necessary for zero range forces. 
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The Skyrme force is usually parameterized as 

V(l,2) = t (l + x P (T )6+lt 1 (l + x 1 P (T )(k' 2 5 + 5k 2 )+t 2 (l + x 2 P (T )k' -5k 

(1) 

+\h P a (l + x 3 P )6 + iW (tn + <r 2 ) • k' x 5k. 
Following the notation employed in ref. ||, particle and pairing density matrices p and p 
are defined as 

p(raq; r'a'q') = ($ | a\ ta , ql a raq | $), (2) 

r'a'q') = (-2a') ($ | | $). (3) 

Owing to the zero range nature of the Skyrme force, the energy functional depends on only 
the local properties of the density matrices. In HFB theory the energy density 7^(r) can 
be expressed as a function of the particle density p(r), the kinetic energy density r(r), the 
pairing density p(r) and the pairing kinetic-energy density t(v). These are defined as 

/'W=E/ , Hi r<T «)> r(r) =V ■Vj2p( ra( l',r , aq) \r'=r , (4) 

aq aq 

P( r ) = P( ra <?; raq) , f(r) = V • p(raq; r'aq) \ r >=r • (5) 

aq aq 

In nuclear matter of isobaric symmetry, Hamiltonian density is given by 



(6) 



^2 3 1 1 

n = — r + (-to + — hp a )p 2 + — {3t 1 + t 2 (5 + 4x 2 )}pr 
2m 8 16 16 

+ |-t (l - x Q ) + ^s(l - x 3 )p°} p 2 + - xt)pf . 
The trial variational function is taken to be of the BCS type, 

i$) = n{^)+^H T / fci Jio), (?) 

q,k 

where k is wave number vector and k =\ k |. Then the densities are given by 

p=- k 2 v 2 (k)dk + — , ~p=--\ k 2 u(kHk)dk, (8) 

9 f^max . „ 9 r- 2 f^ raax , 

r = — / k*v 2 (k)dk + ^k 5 mili , f = J / k A u(k)v(k)dk. (9) 

7T Jfc m i n 07T 7T Jfc m i n 

In the above equations /c min and fc max are defined in terms of a cutoff energy _E C as 

2 2 2m* (p) 2 .2 2m* (p) 

^min — % ^2 ^c> ^max = % H ^2 -^c (10) 
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where m*(p) is the effective mass of nucleon and hp is the Fermi momentum of normal state. 
In order to express the pairing gap as a function of density, a constraint p = po is imposed in 
minimizing the Hamiltonian density (|G|). For this purpose, we utilize the following Routhian, 

H' = H-\p+ C -{p-p,) 2 . (11) 

In the r.h.s., A is the Lagrange multiplier and the third term having a positive constant c 
is included to stabilize the convergence procedure at negative curvature points. If we put 
u(k) = cos9(k) and v(k) = sin9(k), variation with respect to 9(k) yields an equation, 

OH &H 2 

tan2 w = - m m! d ! ; — :■ (12) 

W + ^-X + c(p- P0 ) 



In solving eqs. (§), (0) and fll~2|), a naive iteration procedure does not work when the pairing 
is weak and special manipulations are necessary to obtain the solution. We confirmed that 
the solution actually gives a minimum by giving the Routhian small variations in 9(k) and 
calculating the Routhian. 

The pairing gap depending on k is defined as 

A(k) = -Y,V pp (\k-k' \)u(k')v(k'). (13) 
k' 

The pair-scattering matrix element V pp (q) is of the form, 

V pp {q) = t (l - x ) + -t 3 (l - x 3 )p a + - Xl )q\ (14) 

o I 

where a term depending on the angles is omitted because it vanishes through integration for 
the angles. 

Figure [I] shows the pairing gap A (A; = fcp) as a function of the density characterized with 
the Fermi momentum kp when £' c =10MeV. This value of E c is chosen so as to correspond 
to widely used 2huj configuration space. The results for five parameter sets of Skyrme forces 
(SIII0, SkM*,f§ SGII0, SkP[§ and SkSC4|g) are compared. The term proportional to 
t 3 in eq. ([!]) is treated as a density dependent force. In SIII this term can be dealt with 



as a three-body force. In that case, however, pairing correlation is considerably reduced: 
the maximum pairing gap decreases by a factor of 1/6 and only the normal state solution is 
obtained for density corresponding to kp > 0.7fm _1 . For comparison the result for Gogny 
Dl|| force, a finite range force which reproduces well the pairing properties of nuclei, is taken 
from ref. [|ll||. The figure shows: (i) the gap parameter gives the maximum at kp ~ 0.8 
fm -1 , which is equivalent to about one fifth of the saturation density (kp ~ 1-33 fm -1 ), (ii) 
at the saturation density, with Gogny force the pairing gap of ~ 0.6 MeV still remains while 
with Skyrme forces only the normal state solution is obtained. 

One should notice that, as a result of the zero range nature of the Skyrme force, pairing 
gap may diverge if the summation in eq. QI3D is carried out toward infinity. This implies 
the necessity of a truncation of particle states for the summation. Thus we restricted the 
range of integration in eqs. (§) and Pairing properties of nuclear matter have been 
investigated so far in the context of the neutron star. Jiang and Kuo |L2|| studied pairing in 
nuclear matter by making use of Skyrme forces, but they used the weak coupling approach 
in which pairing gap is constant and not zero only in the vicinity of the Fermi surface. The 
effect of cutoff on the pairing gap has not been studied yet and we extensively discuss it in 
the remainder of this paper. 

Figure |2] shows the pairing gap for kp = 0.8fm _1 as a function of the cutoff parameter 
E c . The pairing gap depends strongly on the cutoff around £' c =10MeV which is employed in 
the calculation shown in Fig. [I]. At larger cutoff, however, there is a plateau for each force. 
The position of the plateau alters from one parameter to another and seems to be associated 
with the momentum transfer q at which the pair-scattering matrix element V pp (q) changes 
the sign. Beyond the plateau the gap starts to increase again. 

We should notice that at high momentum transfer the pair-scattering matrix element 
becomes positive and hence increasing cutoff leads to including a large amount of repulsive 
components in the pairing interaction. Thus increasing cutoff would not be expected to 
increase pairing gap. Nevertheless the calculations show that the pairing gap keeps increasing 
together with the cutoff. It may seem strange at first sight but can be understood as follows. 
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Firstly let us explain the mechanism in which increasing cutoff gives rise to an unlimitedly 
increasing gap. The wi>-factors characterizing the BCS solution are parameterized with 6{k). 
For large cutoff k max , 8(k) takes on negative values for k larger than a certain value (see 
eq. fll2|) ) . In consequence, pairing kinetic-energy density f changes the sign whereas pairing 
density p remains negative, because contribution from high momentum in f is larger than 
that in p. Thus the term proportional to pf in the energy functional (|6]) comes to have the 
effect of lowering energy, notwithstanding the fact that this term usually plays the role of 
raising energy. 

Secondly we consider the physical meaning in such solutions. A similar circumstance, 
in which a repulsive interaction contributes to the pairing correlation, is reported [IT3J] in 
the Relativistic Hartree-Fock-Bogoliubov (RHFB) theory in which the force parameters are 
determined to simulate the Gogny force. They stated that this can happen only in cases 
where the gap parameter as well as the interaction changes the sign as a function of the 
momentum. In the case of the Skyrme force the gap parameter actually changes the sign 
at high momentum and contributes coherently to pairing correlation. But V pp (q) increases 
rapidly as a function of the momentum whereas in the case of finite range force V pp (q) 
diminishes. Therefore the contribution of repulsive part of the force is fictitious in the 
case of the Skyrme force although it is probably realistic for RHFB. The former force is 
constructed so as to describe the ground state properties of nuclei, and corresponds to the 
lowest order expansion in momentum. Its behavior at large momentum transfer does not 
make sense. 

If one takes the cutoff in the middle of the plateau, the dependence of the gap on the 
cutoff can be made reduced considerably. Thus the parameterization (including cutoff) is 
expected to be independent of nuclides. However, since the plateau extends over up to 150 
MeV the question arises as to how to determine the cutoff parameter. Here we propose to 
determine cutoff such that A(A; max ) = 0. 

We also study the dependence of thus-defined cutoff momentum on the density. SIII and 
SkSC4 preserve almost constant value over the whole density whereas the others, especially 
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SkP, possess notable dependence. In the application of HFB to finite nuclei, the pairing gap 
contributes most at the nuclear surface. Therefore we decided to take the cutoff momentum 
to be that at k-p ~ 0.8 fm _1 where the pairing gap has maximum in Fig. [I]. Table 1 shows 
the cutoff momentum and corresponding cutoff energy related through eq. ([TOf ) . 

Figure |3| shows the pairing gap A(k = k F ) calculated with cutoff given in table 1 as a 
function of the density for fcF=0.8fm _1 . Some features are different from Fig. |l|. The peak 
is about three times as high as that in fig. [I] for SkP, SGII and SkM*, whereas that for SIII 
and SkSC4 remains small. The order of the size of gap changes. The position of the peak 
moves to the lower kp (0.6 fm _1 ) for SGII, SkM* and SkP. 

In conclusion, we have investigated the pairing properties of the Skyrme force in nuclear 
matter on the basis of HFB. In particular, the dependence of pairing gap on the cutoff is 
studied. We have observed the appearance of a plateau in the pairing gap at a certain cutoff 
characteristic of each force. The formation of the plateau is attributed to the fact that the 
pair-scattering matrix element changes the sign at a certain transferred momentum. We 
propose to choose the cutoff parameter in the middle of the plateau, so that the pairing gap 
can be made insensitive to the cutoff parameter and the parameterization is more likely to 
be independent of nuclides. 
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TABLE 



TABLE 1. Cutoff momentum & max defined such that A(/c max ) = and cutoff energy E c related 
with k mSuX through eq. ( |10D for fcF=0.8fm _1 . 









£ c [MeV] 


SIII 


1.547 


39.05 


SkM* 


1.857 


61.64 


SGII 


2.076 


80.64 


SkP 


2.966 


169.15 


SkSC4 


1.255 




19.54 
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FIGURE CAPTIONS 



Fig. 1. Density dependence of the pairing gap for Skyrme forces and Gogny Dl force. Here 
the cutoff energy is fixed as E c = lOMeV. 

Fig. 2. The dependence of the pairing gap on the cutoff energy for various versions of Skyrme 
forces. Diamonds indicate the cutoff energy determined such that A(/c max ) = 0. 

Fig. 3. Same as fig. |l]but with the cutoff momentum given in table 1. 
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